Correlated magnetic noise in global networks of gravitational-wave interferometers: 

observations and implications 
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One of the most ambitious goals of gravitational-wave astronomy is to observe the stochastic 
gravitational-wave background. Correlated noise in two or more detectors can introduce a systematic 
error, which limits the sensitivity of stochastic searches. We report on measurements of correlated 
magnetic noise from Schumann resonances at the widely separated LIGO and Virgo detectors. We 
investigate the effect of this noise on a global network of interferometers and derive a constraint 
on the allowable coupling of environmental magnetic fields to test mass motion in gravitational- 
wave detectors. We find that while correlated noise from global electromagnetic fields could be 
safely ignored for initial LIGO stochastic searches, it could severely impact Advanced LIGO and 
third-generation detectors. 
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Introduction. — One of the science goals of 
gravitational-wave (GW) astronomy is to measure 
the stochastic gravitational- wave background (SGWB). 
In cosmological models, the SGWB can be created from 
inflationary physics [1-4], cosmic strings [5, 6], and 
pre-Big-Bang models [7, 8]. Detection of the SGWB 
would constitute a remarkable discovery and could offer 
a unique probe of the history of the universe back to the 
earliest moments after the Big Bang [4]. The SGWB can 
also be created from the superposition of many astro- 
physical objects such as highly magnetized stars [9, 10], 
young or spun- up neutron stars [11-14], core collapse 
supernovac [15], white dwarf binaries [16], super-massive 
black hole binaries [17], and perhaps most promisingly, 
compact binary coalescences (CBCs) [18-20]. Recent 
work [21, 22] suggests that the SGWB from CBCs 
may be detectable with Advanced LIGO (aLIGO) and 
Advanced Virgo [23, 24]. 

The standard procedure for searches for the SGWB 
is to cross-correlate strain data channels from two detec- 
tors [25-29]). By integrating data obtained over a 1 yr- 
long run, it has been possible to achieve astrophysically 
interesting results. The initial LIGO [30] and Virgo [31] 
experiments, which took data from 2000-2010, yielded 
limits on the energy density of GWs that surpass indirect 
bounds from Big Bang Nucleosynthesis and the cosmic 
microwave background [26]. The aLIGO experiment [23], 
scheduled to begin taking data in 2015, is expected to im- 
prove on past results by a factor of « 10 4 [21, 32]. 

A key assumption underlying the SGWB search strat- 
egy is that the noise in each detector is uncorrclated [25, 
33, 34] . This is most easily achieved in spatially separated 
interferometers, and previous analyses [26-29, 35] using 
the LIGO Hanford (LHO) and LIGO Livingston (LLO) 
observatories — separated by a distance of ps 3000 km — 
exhibited no evidence of correlated environmental noise. 
(Self-inflicted correlated noise artifacts, such as a 1 Hz 



comb from similar electronics components at each site, 
have previously been identified and notched.) As GW 
interferometers become more sensitive, however, it is pos- 
sible for subtle global phenomena to produce correlated 
noise at problematic levels. 

In the absence of correlated noise, the sensitivity of a 
search for the SGWB is limited only by the observation 

1/2 

time t bs, with the signal-to- noise ratio growing like i {, 8 . 
Correlated noise, however, produces a systematic error, 
which cannot be reduced through integration. To the 
extent that it cannot be mitigated through instrumental 
(re)design and/or background subtraction, it constitutes 
a fundamental limit for SGWB searches. 

Previous work [25, 33, 34] has identified Schumann 
resonances as a potential source of correlated noise 
for widely separated detectors, especially for second- 
generation experiments such as aLIGO. Schumann res- 
onances, predicted in 1952 [36, 37] and observed soon 
thereafter [38, 39], are global electromagnetic (EM) res- 
onances in the cavity formed by the surface of the Earth 
and the ionosphere. The cavity is excited by a back- 
ground of w 100 lightning strikes per second around the 
world with 20 — 30 kA of current and lengths of 3 — 5 km. 
The resonances produce magnetic fields on the Earth's 
surface of 0.5 - l.OpTHz" 1 / 2 [40]. Observed in the 
time domain, 10 pT bursts appear above a 1 pT back- 
ground [41] at a rate of pa 0.5 Hz. The primary reso- 
nance is at 8 Hz with secondary and tertiary harmonics 
at 14 Hz and 20 Hz respectively. The peaks exhibit a 
spectral width of sa 20% and vary seasonally and with 
proximity to lightning storms. 

Global magnetic fields such as Schumann resonances 
can cause correlated noise in gravitational- wave interfer- 
ometers by inducing forces on magnets mounted directly 
on the test masses for position control (as in initial LIGO) 
or on magnets (and magnetically susceptible material) 
higher up the test-mass suspension system (as in aLIGO). 
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Since Schumann resonances produce coherent EM fields 
over sa 1000 km distances, it may be possible for them to 
introduce correlated strains in widely separated detec- 
tors. Recent observations of Schumann resonances and 
a discussion of their properties are given in [40, 42, 43]. 
Other EM phenomena such as solar storms, currents in 
the van Allen belt [44] , and anthropogenic emission may 
also contribute to correlated EM noise. 

In this Letter we report on measurements of correlated 
EM noise (including contributions from Schumann reso- 
nances) observed in magnetometers situated at widely 
separated GW observatories. Combining these results 
with prior measurements of the coupling between mag- 
netic fields and LIGO test mass motion, we infer the level 
of correlated strain noise present in initial LIGO SGWB 
analyses. Then we derive magnetic isolation specifica- 
tions for the next generation of advanced detectors to 
ensure that future SGWB analyses are not limited by cor- 
related noise. We conclude by discussing the implications 
for third-generation detectors and possible strategies for 
subtracting correlated noise. 

Formalism. — Searches for the SGWB typically seek to 
measure the logarithmic energy density of GWs 
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where / is GW frequency, p c is the critical energy density 
required for a closed universe, and dpew is the energy 
density of GWs between / and f + df. An optimal es- 
timator for f^GW, integrated over some detection band, 
can be constructed from the strain time series of two GW 
interferometers [25]: 



Y(f) = —R C (Q(m(f)S 2 (f)) 



(2) 



Here Q(f) is a filter function (see [25]), ST is the data 
segment duration, and s/(/) is the Fourier transform of 
the strain series measured by detector /. By combining 
data from many data segments and frequency bins, it is 
possible to construct an optimal broadband estimator for 
the entire science run, (see [25] for details). 

Following [25], the strain observed by each detector 
contains contributions from signal and noise and can be 
written as 



si(f) = h(f)+ni(f), 



(3) 



where hi(f) is the GW strain induced in detector / and 
n/(/) is the noise. If n\(f) and ?t-2(/) are uncorrelated, 
an optimal filter Q(f) can be chosen such that 



(Y(f)) = n GW (/), 



(4) 



where the angled brackets denote an expectation value. 
An estimator for the uncertainty of Y(f), denoted ay(f), 

— 1/2 

grows like i obs . In the absence of correlated noise, Y(f) 



is an unbiased estimator and the sensitivity of the search 
is limited by i bs- 

In the presence of correlated noise, Eq. 4 becomes 



where 



(Y(f)) = ficwCf) + fM/) 



Ojvee— Re[Q(/)(^(/)n 2 (/))] 



(5) 



(6) 



is the correlated noise. represents a system- 

atic bias. If f2jv(/) ^ ®Y{f), then the search is lim- 
ited by how well f2jv(/) can be estimated. If, however, 
^at(/) <C then the correlated noise term may be 

safely ignored. 

Before we continue, we define coherence, which is use- 
ful for determining if two channels are correlated. The 
coherence between channels 1 and 2 is given by: 



coh(/) 



n(f)Uf) 
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(7) 



Here the overline denotes time-averaging over N seg- 
ments. If Si and S2 are independent, Gaussian, and sta- 
tionary random variables, then (coh(/)) = 1/N. Devia- 
tions from coh(/) ps 1/JV are evidence that one of these 
three assumptions is violated. Integration over long ob- 
servation periods probes below the uncorrelated noise to 
unearth low-level coherent features. 

It is also useful to define coherence in terms of cross- 
and autopower 
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where 



is the autopower for detector I and 



812(f) = jf§l(f)h(f) 



(8) 



(9) 



(10) 



is the cross-power for the detector pair 1 and 2. M is 
a discrete Fourier transform normalization constant. We 
respectively refer to y/ Si(f) and y/\Si2(f)\ as the auto- 
and cross- amplitude spectra. 

Procedure and results. — The LIGO and Virgo obser- 
vatories are equipped with sensors to measure environ- 
mental noise that may affect strain measurements [30]. 
Of particular interest here are LIGO's Bartington Mag- 
03 magnetometers and Virgo's Metronix MFS-06 magne- 
tometers, both housed inside observatory buildings. We 
calculate coh(/) for magnetometer channels at different 
detectors to determine if they exhibit excess coherence. 
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Data are broken into 10 s segments yielding 0.1 Hz resolu- 
tion. Example multi-month coherence spectra are shown 
in Fig. 1 for LHO-LLO and LHO- Virgo magnetometer 
pairs. There are several noteworthy features. 

First, we observe Schumann peaks (marked with black 
circles) , three of which are coherent between all three ob- 
servatories. The peak frequencies / = 7.8, 13.9, 20.5 Hz 
etc. are consistent with expected values. While Fig. 1 
shows just two spectra, we observe similar peaks for 
channel pairs associated with different magnetic field di- 
rections [49], for magnetometers located at different lo- 
cations at the observatories, and for data taken during 
different observing periods. By plotting the cumulative 
coherence as a function of time, we determine that the 
observed peaks are not due to non-stationarity effects (a 
small number of data segments with anomalously high 
power) . We observe no broadband coherence at high fre- 
quencies / > 200 Hz for LHO-LLO and / > 50 Hz for 
LIGO- Virgo. 

Second, we observe narrowband features (marked with 
green diamonds) at frequencies that are expected for 
noise from electronics systems. The LHO-LLO coher- 
ence spectrum exhibits an expected line at 60 Hz from 
the electrical power systems. The LHO- Virgo spectrum 
exhibits lines at 16 Hz, 50 Hz (extending beyond the plot 
range), and 100 Hz. The lines at 50 Hz and 100 Hz are 
produced by nearly monochromatic electrical power sys- 
tem frequencies at Virgo aligning in frequency with mem- 
bers of a 10 Hz comb produced by a timing synchroniza- 
tion system at LHO. Since we are concerned here with 
magnetic fields that are correlated over global distance 
scales, we subsequently remove these line artifacts. 

Finally, as an additional cross-check, we repeated the 
coherence calculation while shifting one time series 10 s 
with respect to the other. The time-shifted spectra ex- 
hibit no broadband excess while electronic lines remain, 
which supports our conclusion that the broadband peaks 
are due to global magnetic correlations. 

Implications. — In order to investigate the effect of EM 
correlated noise during initial LIGO, we utilize measure- 
ments of the coupling function T(f) (the absolute value 
of the transfer function) , which describes the induced test 
mass motion per unit environmental magnetic field. The 
measurements, performed multiple times during each sci- 
ence run, utilize magnetic injection coils (placed outside 
but nearby the vacuum chambers) to create alternat- 
ing magnetic fields. The injection coils are located at 
a distance great enough that the magnetic field gradi- 
ents at the test mass are dominated by the effects of 
local field-altering structures. The resultant test mass 
motion is measured using the GW strain channel. The 
initial LIGO magnetic coupling functions for each test 
mass at both sites were found to be within a factor 
of two of T(f) = 0.003(./yiHz)- 3 m/T [45]. Also, the 
coupling function was found to be the same (to within 
a factor of two) for magnetic fields injected both per- 
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FIG. 1: Magnetometer coherence spectra for LHO-LLO dur- 
ing the LIGO S5 science run (top, t b s = 330 dy) and for LHO- 
Virgo during S6-VSR2/3 (bottom, t b s = 100 dy). Schumann 
resonance peaks are indicated with black circles while elec- 
tronic noise lines are indicated with green diamonds. The red 
dashed line indicates the average value expected for uncor- 
rected noise. Some LHO-LLO peaks are obscured by 60 Hz 
electronic noise. The frequency resolution is 0.1 Hz. 

pendicular and parallel to the beam axis because of lo- 
cal field-altering structures [45]. We therefore proceed 
with the simplifying assumption that the coupling func- 
tion is approximately the same at each interferometer: 
Ti(/) = T 2 [f) = T(f). Further, in the analysis that 
follows, we conservatively use the magnetometer spec- 
tra exhibiting the strongest coherence at each LIGO site. 
The least coherent spectra are approximately an order of 
magnitude less coherent, which would yield weaker cou- 
pling function requirements by ~ 2. 

The strain noise cross-power spectrum induced from 
magnetic fields is 

JVS(/) « ^(/)T 2 (/) |M 12 (/)| , (11) 

where 



M 12 (/) = -mi(f)m 2 (f) (12) 

is the magnetometer cross-power spectrum (see Fig. 2), 
rhi(f) is data from magnetometer channel /, and L is 
the interferometer arm length. By relating Nf^(f) to the 
absolute value of Mi 2 (/), we conservatively assume that 
the phase of M i2 (/) yields the worst possible contamina- 
tion in the strain channel. The superscript m in Eq. 11 
labels this noise spectrum as correlated strain noise of 
magnetic origin as opposed to uncorrelated strain noise. 

In the absence of correlated noise, the sensitivity of a 
stochastic search is determined by the uncorrelated noise 
spectrum, estimated by a weighted average of the product 
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FIG. 2: Magnetic cross-amplitude spectra \/|Mia(/)| during 
the S5 and S6-VSR2/3 LIGO- Virgo science runs. Note that 
Fig. 2 and Fig. 1 are related through Eq. 8. The frequency 
resolution is 0.1 Hz. The blue S5 HL A and violet S6 HL 
traces utilize the same magnetometer pair, whereas the black 
S5 HL B trace comes from a magnetometer pair only available 
during S5. 



of the autopower spectra [25]: 



N? 2 (f) 



-1/2 



(13) 



The superscript u in Eq. 13 labels this noise spectrum as 
uncorrelated noise, which decreases with t ^ s . 

In Fig. 3 we show the correlated strain noise ampli- 
tude spectrum y/N$(f) for the initial LIGO 4 km HI 
and LI detectors located at LHO and LLO respectively 
(red). Alongside, we plot \/iV" 2 (/) for aLIGO at design 
sensitivity with 1 yr of coincident data (black), and for 
initial LIGO with 330 dy of coincident data (purple). We 
estimate that the correlated noise was below the uncorre- 
lated noise during initial LIGO, which is consistent with 
back-of-the-envelope estimates from [25], previous noise 
studies [46] , and past observational results [26] . 

The aLIGO experiment, however, is expected to 
achieve a factor of 10 better strain sensitivity than ini- 
tial LIGO with a wider detection band going as low as 
10 Hz. Using Eq. 11, and calculating Nf 2 (f) for aLIGO 
design sensitivity, we can solve for T(f) — the coupling 
function at which the correlated noise becomes compa- 
rable to the uncorrelated noise. T(f) is shown in Fig. 4 
for two different observation durations. We also include 
T(f) calculated for the LHO- Virgo detector pair. 

The coupling function in Fig. 4 may not be conser- 
vative enough because even when the correlated noise 
is below the uncorrelated noise, it can become signifi- 
cant through integration if it occurs in many frequency 
bins. The enhancement of sub-threshold signals grows 
like the square root of the number of frequency bins. 
Given a 100 Hz band of correlated noise and a resolu- 
tion of 0.25 Hz, it would be prudent to aim for a factor 



FIG. 3: Strain amplitude spectra for correlated and uncor- 
related noise. Black is the uncorrelated noise -J 7V" 2 (/) for 
the H1L1 detector pair operating at Advanced LIGO design 
sensitivity and assuming 1 yr of integration. Purple indicates 
the uncorrelated noise \JN^ 2 (f) achieved during initial LIGO 
using w 300 dy of coincident data. Red is */ (/) (the esti- 
mated correlated noise due to EM fields during initial LIGO). 
Electronic noise lines have been notched. The spectra have 
been scaled to assume a frequency resolution of 0.25 Hz, which 
is typical for stochastic searches [26]. 




FIG. 4: The design-sensitivity Advanced Virgo and aLIGO 
critical coupling function T(f) above which correlated noise 
is comparable to uncorrelated noise. The critical coupling 
function is calculated by considering the sensitivity obtained 
from a given pair of detectors, each assumed to have the same 
coupling function. Electronic noise lines have been notched. 
We also include the critical coupling function (solid black) for 
initial LIGO. The measured initial LIGO coupling function 
(dashed black) falls below it. 



of v400 ss 20 lower than the coupling shown in Fig. 4. 

Correlated magnetic events may also be of concern 
for GW burst searches. A 10 pT magnetic event and 
a coupling of 1 x 10 _7 m/T would generate a strain 
of h w 2.5 x 10 -22 . However, the limited coherence 
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coh(/) 10 -3 would likely preclude coincident false de- 
tections. Magnetic bursts and their influence on the in- 
terferometers are the subject of ongoing research. 

The best means of guarding against EM-induced cor- 
related noise is to minimize magnetic coupling to the 
test masses, e.g., through the removal of magnetic com- 
ponents. As future interferometers, e.g., [47], look to 
further reduce (uncorrelated) strain noise and to oper- 
ate at lower frequencies, the required coupling function 
constraints will become increasingly stringent, and other 
remedies may be necessary. 

One possibility worthy of study is to use magnetome- 
ters to subtract the correlated noise, e.g., with a Wiener 
filter scheme. Similar tactics are under investigation to 
subtract Newtonian gravitational noise in GW interfer- 
ometers [48] . Effective subtraction requires precise mea- 
surement of global EM fields, which can lurk underneath 
an order of magnitude stronger local fields. Future stud- 
ies are necessary to determine the residual contamination 
expected from realistic background subtraction. 
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